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ABSTRACT 

An ionizing UV background (UVB) inhibits gas accretion and photo-evaporates gas from the 
shallow potential wells of small, dwarf galaxies. During cosmological reionization, this effect 
can result in negative feedback: suppressing star-formation inside HII regions, thus imped- 
ing their continued growth. It is difficult to model this process, given the enormous range of 
scales involved. We tackle this problem using a tiered approach: combining parameterized 
results from single-halo collapse simulations with large-scale models of reionization. In the 
resulting reionization models, the ionizing emissivity of galaxies depends on the local values 
of the reionization redshift and the UVB intensity. We present a physically-motivated ana- 
lytic expression for the average minimum mass of star-forming galaxies, M m ; n , which can be 
readily used in modeling galaxy formation. We find that UVB feedback: (i) delays the end 
stages of reionization by Az < 0.5; (ii) results in a more uniform distribution of HII regions, 
peaked on smaller-scales (with large-scale ionization power suppressed by tens of percent); 
and (iii) suppresses the global photoionization rate per baryon by a factor of < 2 towards the 
end of reionization. However, the impact is modest, since the hydrodynamic response of the 
gas to the UVB occurs on a time-scale comparable to reionization. In particular, the popular 
approach of modeling UVB feedback with an instantaneous transition in M m ; n , dramatically 
overestimates its importance. UVB feedback does not significantly affect reionization unless: 
(i) molecularly-cooled galaxies contribute significantly to reionization; or (ii) internal feed- 
back processes strongly couple with UVB feedback in the early Universe. Since both are 
considered unlikely, we conclude that there is no significant self -regulation of reionization by 
UVB feedback. 

Key words: cosmology: theory - early Universe - galaxies: formation - high-redshift - 
evolution 



1 INTRODUCTION 

Early generations of luminous objects reionized the Universe 
within the first ~billion years following the Big Bang. This pro- 
cess is expected to be fairly extended and inhomogeneous. As ion- 
izing radiation spread out from early galaxies, it heated the inter- 
galactic medium (IGM), affecting its cooling properties, and photo- 
evaporated gas from shallow potential wells. In this picture, galax- 
ies forming inside an already-ionized IGM would have a depleted 
gas reservoir (and by extension, fewer stars) compared with those 
forming inside the neutral IGM. Therefore, the ionizing ultraviolet 
background (UVB) during reionization results in a negative feed- 
back mechanism, hindering the further growth of ionized regions. 
This negative feedback could delay reionization, as well as result in 
a more uniform reionization morphology. Therefore understanding 
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UVB radiative feedback is crucial in developing models of reion- 
ization, as well as interpreting observations. 

The main difficulty in self-consistently modeling UVB feedback 
lies in the enormous dynamical range required. Ideally, we would 
need to fully resolve the internal structure of the dominant popu- 
lation of atomically-cooled halos (with virial radii ~ 1 proper kpc; 
IBarkanalsT Loeb 2001 ) in cosmological simulations which resolve 
ionizing structure on > 100 comoving Mpc scales (e.g. Furlanetto 
|et al.|[2004| . Even state-of-the-art single galaxy simulations must 
employ some analytic prescriptions in estimating internal feed- 
back (e.g. Hopkins et al. 201 I ): therefore compromises must be 
made. Generally one of two approaches are followed. The first ap- 
proach uses relatively small simulation boxes (~ few cMpc on a 
side), attempting to model the sub-grid interstellar medium (ISM) 



physics, through analytical prescriptions (e.g. 


Petkova & Springel 


2011||Finlator et al.|2011||Wise et al.|2012||. 


ohnson et al. 2012; 


Hasegawa & Semelin 2012). Here it is important to explore reason- 



able extremes of all sub-grid recipes, and ensure resolution con- 
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vergence. Even so, it is currently impossible to model the relevant 
scales of reionization. The second approach simulates reionization 
on large-scales, but does not model the internal structure of galax- 
ies; instead, simple analytic (and relatively ad-hoc) prescriptions 
are used to account for UVB feedback (e.g. Iliev et al. 2007; Mc- 
IQuinn et al.|2007l l. 

In this paper we follow the later approach. However, rather than 
adopting ad-hoc prescriptions, we base our treatment of UVB feed- 
back on results from parameter studies of ID collapse simulations 
(Sobacchi & Mesinger 2013 hereafter Paper I). In Paper I, we ob- 
tained expressions for the baryon content of galaxies, residing in 
regions which were reionized at redshift zm. In this work, we in- 
clude these expressions in semi-numerical simulations of inhomo- 
geneous reionization. We present large-scale reionization simula- 
tions which include a physically-motivated prescription for UVB 
feedback. With these simulations, we study the importance of UVB 
feedback in regulating reionization. 

This paper is organized as follows. In [j2]we describe our UVB 
feedback prescription, and our semi-numerical reionization simula- 
tions. In *j3]we present our results, comparing different inhomoge- 
neous reionization models. In !|4]we discuss our results and present 
our conclusions. Throughout we assume a flat ACDM cosmology 
with parameters (fi m , fi A , Mb, h, <t 8 , n) = (0.27, 0.73, 0.046, 0.7, 
0.82, 0.96), consistent with WMAP results ( |Komatsu et al.|2011| >. 
Unless stated otherwise, we quote all quantities in comoving units. 



redshift-dependent halo mass (e.g. |Barkana & Loeb|200T) : 
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where /i is the mean molecular weight, Sl^ — 
O m (1 + zf I [fi m (1 + zf + Oa] and A c = 187r 2 +82d-39d 2 
with d = fi£ - 1. 

Fitting to a broad parameter space of ID collapse simulations, 
in Paper I we presented the following expression for M cr it (techni- 
cally defined as the total halo mass at which the baryon fraction is 
1/2 of the global value, Sit /Sim)'- 
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with best-fit parameters: 

(M ,a,b,c,d) = (2.8 x 10 9 M Q , 0.17, -2.1, 2.0, 2.5) (4) 

The critical mass depends on (i) redshift z; (H) the redshift zm 
when the halo is exposed to a UVB; (Hi) the UVB intensity J21 Q 

J (v) = J21 iylvYL)'" x 10~ 21 erg s" 1 Hz -1 



(5) 



where vh is the Lyman limit frequency and a = 5 corresponds to 
a stellar-driven UV spectrum (e.g. |Thoul & Weinberg 1996). 



2 MODELING REIONIZATION WITH UVB FEEDBACK 
2.1 Minimum halo mass to host galaxies 

In the simplest UVB-regulated scenario^ there are two fundamen- 
tal halo mass scales regulating the gas reservoir available for star 
formation: (i) M cr it, the characteristic or critical mass below which 
baryons are photo-evaporated or cannot efficiently accrete onto 
their host halos; (ii) M coo i. the cooling threshold required for gas to 
efficiently cool, collapse and form stars inside the halo. Therefore, 
we can define the minimum mass of star-forming halos as: 

M min = max [Mcooi, Merit] • (1) 

Here we mainly consider the advanced stages of reionization, 
where cooling through molecular hydrogen is expected to be rel- 
atively inefficient due to a disassociating background (e.g. |Haiman| 
|et al.|[T997) l. This corresponds to an atomic cooling threshold of 
a constant virial temperature, T V1I ~ 10 4 K, or analogously a 



1 We caution that we do not directly include internal feedback processes, 
which, depending on the ISM treatment, could be important in determin- 
ing the star formation rate (SFR) of dwarf galaxies at high-redshifts (e.g. 
|Pawlik & Schaye|2009| |Finlator et al.|2011| |Wyithe & Loeb|2012) . Al- 
though winds can push gas to the outskirts of galaxies, where it can be 
more easily photo-evaporated by a UVB, this amplification effect is mod- 
est at high-redshifts (at redshifts not far below 2jn; e.g. Pawlik & Schaye| 
2009 , Finlator et al. 201 1 1. Hence, during most of reionization, it is probably 
reasonable to model internal feedback independently from UVB feedback 
(e.g. by decreasing the ionizing efficiency parameter, £, below). Since we 
are primarily concerned with the relative impact of UVB feedback, we do 
not expect internal feedback to alter our main conclusions. 



2.2 Semi-numerical reionization simulations 

To model cosmological reionization, we use a parallelized version 
of the publicly available semi-numerical simulation, 21CMFASlj^] 
We generate the IGM density and source fields by: (i) creating a 
3D Monte Carlo realization of the linear density field in a box with 
sides L = 300 Mpc and N = 1600 3 grid cells; (ii) evolving the 
density field using the Zel'dovich approximation (Zeldovich 19701, 
and smoothing onto a lower-resolution N = 400 a grid; (Hi) using 
excursion-set theory jPress & Schechterj|1974| |Bond et aL||1991| 
|Lacey & Cole|1993||Sheth & Tormen|1999) on the evolved den- 
sity field to compute the fraction of matter collapsed in halos big- 
ger than M m i n , thus contributing to reionization (see Mesi nger &| 
|Furlane tto 2007 ; Mesinger et al.|2011| for a more detailed descrip- 
tion of the code). 

The ionization field is computed by comparing the integrated 
number of ionizing photons to the number of baryons, in spheri- 
cal regions of decreasing radius R (i.e. following the excursion-set 
approach of Furlanetto et al. 2004). Specifically, a cell located at 
spatial position and redshift, (x, z), is flagged as ionized if: 

£/coii(x,i?,z,M min ) > 1 (6) 

where / co n (x, ii, z, M m i n ) is the fraction of collapsed matter in- 
side a sphere of radius R residing in halos larger than M m i n , and 

2 In Paper I, we ignored self-shielding when computing the critical mass. 
The importance of self-shielding on such dwarf galaxies at high-redshifts is 
not well known, with some authors suggesting that it does not have a large 
impact on the relevant photo-evaporation time-scales (e.g. Shapiro et al. 
|2004| |Uiev et al.|2005) . Nevertheless, it should be noted that neglecting 
self-shielding overestimates the impact of UVB feedback on reionization, 
making our main conclusions conservative. 

3 http://homepage.sns.it/mesinger/Sim 
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Figure 1. Slices through our ionization boxes, with HII regions in black and HI regions in white. From top to bottom we show the runs: NF, IF, J=0.01, SC, 
J=l. 



£ is an ionizing efficiency, defined below. Starting from J? m ax = 
■Rmfp = 30 Mpc (roughly corresponding to the ionizing photon 
mean free path in the ionized IGM, i? m f p , at z ~ 6 So ngaila"&| 
|Cowie 2010; McQuin n et al.|20"TT] > the smoothing scale, 7?, is de- 
creased, and the criterion in eq. ([6]l is re-evaluated. At the cell size, 
the partial ionizations from sub-grid sources are evaluated, includ- 
ing Poisson noise around the mean value of the collapse fraction 
(e.g. |Mesinger et al.||2011) . This algorithm results in ionization 
fields which are in good agreement with cosmological radiative 
transfer algorithms ( Zahn et al. 201 1 1. 

The ionizing efficiency can be written out as: 



i = 30 
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where / csc is the fraction of UV ionizing photons that escape into 



the IGM, N-, is the number of ionizing photons per stellar baryon, 
/* is the fraction of galactic gas in stars, and n TCC is the mean num- 
ber of recombinations per baryon in the IGM. Our fiducial choice 
of £ = 30 results in an ionization history in agreement with the 
WMAP observed value of r e {Komatsu et al. 2011); however we 
also explore other values below. Although our models depend only 
on the product in eq. (|7|, we show on the RHS some reasonable 
values for the component terms. iV 7 « 4000 is expected for PopII 
stars (e.g. Barkana & Loeb|2 005 1, and studies of the high-redshift 
Lyman alpha forest suggest n re c ~ in the diffuse IGM (Miralda- 
|Escude|2003||Bolton & Haehnelt|20071|McQuinn et al.|201 1| note 
that we instead incorporate photon sinks in the form of a mean 
free path through the ionized IGM, governed by the separation 
of Lyman limit systems, LLSs). On the other hand, the parame- 
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ters /esc and /» are extremely uncertain in high-redshift galaxies 
(e.g. |Gnedin et al.|2008[|Wise & Cen|2009l|Ferrara & Loeb|2012| >, 
though our fiducial choices are in agreement with high-redshift 
galaxy luminosity functions (e.g. Robertson et al. 2013}. 

The main improvement in this work is that we use the spatially- 
dependent value ofMmi n , computed according to eq.^(rather than 
a homogeneous value M m in = M coo i as is commonly done). We 
keep track of the redshift when each cell was first ionized, jzin(x), 
updating the local value of M CI i t (z, zin, J21) as reionization pro- 
gresses (for details on how we compute the inhomogeneous J21, 
see the following section). Neutral cells are assigned the value 
M coo i(z). We then use the average value, (A/ m i n )_R in the ioniza- 
tion criterium of eq. JSJ. The drawback of this addition is that it 
requires "running-down" the simulation from high redshifts, rather 
than being able to independently generate the reionization mor- 
phology at a given redshift. 

2.2.1 Incorporating an inhomogeneous, self-consistent UVB 

We also compute the average intensity inside each HII region, for 
use in eq. QJ. The mean emissivity (number of ionizing photons 
emitted into the IGM per second per comoving volume) can be 
written as: 

r r AT" - /coll ^_ /coll /ox 

e « /./cscA^ns,— — w t,n b —— , (8) 

where n b is the mean baryon number density inside the HII region, 
and t* is the star-formation time-scale (which we take to be t* = 
0.3 x tn). The intensity is proportional to the mean free path inside 
the HII region (taken to be « R) multiplied by the emissivity. If the 
emissivity is spectrally distributed as v~ a (eq.[5}, we can write the 
average UVB intensity in a given HII region of characteristic size 
R (in erg s _1 Hz -1 sr _1 proper cm -2 ) as 

(J 2 i)hii = {1 + z)2 Rhan b ^ . (9) 

Due to the clustering of dark matter halos, the relevant intensity 
at galaxy locations will be higher than the average UVB intensity in 
eq. Therefore when calculating A/ m j n , we use (J2i)haio,HII — 
/jbias x (J2i}hii, with /jbias = 2, consistent with results from 
Mesinger & Dijkstra ( 2008 1. We also note that the simplification of 
a uniform UVB inside each HII region seems to be reasonable (see 
Fig.|7]and the associated discussion). 

Furthermore, this formalism is not fully self-consistent, since 
our expression for M cr it was computed for a time-independent 
UVB. Nevertheless, we find that the mean UVB is roughly constant 
in HII regions, somewhat validating this approximation. More im- 
portantly, feedback is largely insensitive to the exact value of J21 
since the gas is always heated to T ~ 10 4 K (specifically, in Paper 
I we find Mont oc J 017 ). 

2.3 Runs 

The efficient, semi-numerical approach outlined above allows us 
to explore several scenarios of UVB feedback during reioniza- 
tion. As already mentioned, our fiducial model assumes £ = 30, 
^cool = 10 4 K and 7? mfp = 30 cMpc, though we explore other 
values throughout. We have five main prescriptions for implement- 
ing UVB feedback: 

• No feedback (NF): we assume M m in = M coo i, correspond- 
ing to a constant virial temperature T CO oi. This model completely 
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Figure 2. Evolution of the global neutral fraction Shi with different mod- 
els of radiative feedback. Dot-Dashed: no feedback (NF). Short-dashed: 
J21 = 0.01 inside HII regions. Dotted: J21 = 1 inside HII regions. 
Dashed: instantaneous feedback (IF). Long-dashed: self-consistent model 
(SC). 

neglects UVB feedback, as is commonly done in reionization liter- 
ature. 

• Instantaneous feedback (IF): in this extreme model, feed- 
back leads to an instantaneous transition in Mmin as soon as the 
halo is exposed to a UVB, thus ignoring the hydrodynamical re- 
sponse time-scale of the gas. This is the most common prescrip- 
tiorj 4 ] for UVB feedback (e.g. |Wyithe & Loeb||2003[ |Onken &" 
Miralda-Es cude|2004[ [Biev et al.|2007||Wyifhe & Morales|2007[ 
Wyithe & Cen|2007||Kulkarni & Choudhury|2011||Alvarez et al'| 
2012| |Ahn et al.||2012| >, and we take A/ m i n = M coo i in HI re- 
gions and M m i n = 10 9 Mq in HII regions, consistent with pre- 
vious works. The IF and NF runs therefore bracket the expected 
impact of UVB feedback. 

• J=l: we calculate M m i n with eq. [T| using a constant UVB 
intensity J21 = 1 in HII regions. This intensity is larger than esti- 
mates from the Lyman alpha forest at z ~ 6, J21 ~ 0.1 (jBolton 
|& Haehnelt|2007l|Wyithe & Bolton|20"TT] |Calverley et al.|2011| > 
Therefore this model places an upper limit to the impact of UVB 
feedback (which is more physically relevant than the upper limit 
provided by the IF model). 

• J=0.01: like the previous model but with J21 = 0.01. 

• Self -consistent (SC): we calculate Mmin using eq.[T| and we 
use eq.|9]to calculate the inhomogeneous UVB in HII regions. 



3 RESULTS 

3.1 Reionization History 

In Fig. [T| we show slices through the ionization boxes of our five 
fiducial models: NF, IF, J=0.01, SC, J=l (top to bottom). The 
corresponding electron-scattering optical depths are: r c = 0.088, 
t c = 0.073, r e = 0.087, r e = 0.085, r e = 0.08^] consistent at 
la with WMAP results ( |Komatsu et al.|2011) . 

4 In simulations, this generally amounts to using the simulation time-step 
(generally ailO Myr) as the feedback time-scale. 

5 UV-driven reionization is "inside-out" on large-scales, with overdensities 
hosting the biased galaxies being the first to ionize. Mesing er et al.| j2012) 
noted that this correlation between the density and ionization fields results 
in a ft: 4% higher t c than what would be expected assuming the two fields 
are uncorrelated: {xi xn b ) ^ (xi) X (rib). F° r example our self-consistent 
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As expected, UVB feedback delays reionization and results in 
a more uniform distribution of HII regions (e.g. |McQuinn et al.| 
[2007] |Iliev et al.|2007| >. However, our self-consistent model is not 
remarkably different from the one which ignored UVB feedback al- 
together (NF). Most striking is the effect of an unphysical instanta- 
neous transition in M m i n : the IF model dramatically over-predicts 
the importance of UVB feedback even with respect to the extreme 
J=l model. We can therefore immediately surmise that the delay- 
in the hydrodynamic response of the baryons to the UVB is funda- 
mental in properly modeling UVB feedback. 

This qualitative result is confirmed by the evolution of the 
global neutral fraction Shi, shown in Fig. [2] As expected, the self- 
consistent model shows intermediate features between the J=0.01 
and J=l models. Since the transition of M m i n is not instantaneous, 
only the late stages of reionization are delayed. Comparing the SC 
and NF curves, we see that UVB feedback delays the later stages 
of reionization by Az < 0.5. On the other hand, this delay is sig- 
nificantly overestimated by a model with instantaneous feedback: 
Az « 1.5. 



3.2 Reionization Morphology 

We now proceed to quantify the impact of UVB feedback on 
the morphology of reionization. In Fig. [3] we show the ioniza- 
tion power spectrum (A xx = k 3 /(2n 2 V) (|5 xx | 2 )fc, with <5 XX = 
£hi/Shi — 1), while in Fig.|4]we show the size distributions of the 
HII regions]^] Panels correspond to different stages of reionization: 
Shi = 0.75,0.50,0.25, (left to right). We compare our fiducial 
models of UVB feedback: no feedback (yellow dot-dashed line); 
J21 = 0.01 (magenta short-dashed line); J21 = 1 (blue dotted 
line); instantaneous feedback (green dashed line); self-consistent 
(blue long-dashed line). 

As expected from Fig. [T] we see that ignoring feedback over- 
predicts the large-scale ionization power by tens of percent. The 
same trend can be seen in the HII size distributions, with the peak 
in the NF model being shifted by w 50% with respect to the SC 
model. The IF model on the other hand over-estimates these feed- 
back effects by a comparable amount. We therefore confirm earlier 
predictions that UVB feedback results in smaller, more uniform 
ionized regions (Mc Quinn et al.|2 007 1. This is easy to understand: 
the largest HII regions correspond to the most biased locations of 
the density field, and are therefore ionized earliest. It is in these 
regions that UVB feedback has had the most time to quench star- 
formation. Smaller, late-forming HII regions are less affected, and 
therefore more abundant at fixed Shi. 

Nevertheless, the impact of UVB feedback is relatively minor, 
when compared at fixed xm'- the power spectra of all models agree 
to within a factor of < 2. This is due to the fact that the halo 
bias evolves only weakly over the mass scales relevant to UVB 
feedback; hence morphology at a given Shi is relatively robust 



model has r c = 0.085 instead of r c = 0.083 as would be expected ignor- 
ing this correlation. 

6 We calculate the size distributions following the prescription in |Mesinger| 
& Furlanetto 1 2007 1. Namely, we randomly choose an ionized cell and then 
record the comoving distance to the edge of the HII region in a random 
direction. Then we define the size distribution as Rdp/dR, where p is the 
probability that this distance is between R and R + dR. Technically we 
define the edge of an HII region as the first cell with :thi > 0.5. This 
choice of threshold can affect the size distributions somewhat; however, we 
are interested in the relative difference between the models, which is much 
more robust to this choice (Friedric h et al.|201 1| . 



(McQ uinn et al.||2007| >. In other words, reionization morphology 
is much more sensitive to Shi than it is to UVB feedback effects. 
This implies that previous predictions for reionization morphology 
(at fixed Shi) are not particularly sensitive to this source of astro- 
physical uncertainty. 

3.3 Evolution of M m i„ 

The results of the previous sections can be more readily understood 
by looking at the redshift evolution of the average value of the min- 
imum halo mass hosting star-forming galaxies, M m i n , shown in 
Fig. [5] using the same line styles as above. This quantity is also 
fundamental in semi-analytical models of galaxy formation. 

By construction, the M m i n curve in the NF model simply fol- 
lows the Tvh- = 10 4 K isotherm. The curve corresponding to 
the SC model is closer to the J=l model than the J=0.01 one 
(see Fig. [7] and associated discussion). It begins "pealing off" the 
T v i r = 10 4 K isotherm at z ~ 10, when reionization is already 
well underway (Shi ~ 0.6). In contrast, M m i n in the IF model 
immediately starts increasing as reionization progresses, since in 
that model the average minimum halo mass is simply M m i n = 
Sm-Mcooi + (1 — Shi)1O 9 A/0. Because it ignores the delay asso- 
ciated with UVB feedback, the IF model over-predicts M m i n by a 
factor of ~2— 3, with respect to the SC model. 

For comparison we also show the evolution of M m i n assuming 
a homogeneous reionization at zm = 9.3 = z rc , corresponding 
to the redshift when Shi = 0.5 in the self-consistent model (ma- 
genta solid line). This formula under-predicts M m i n by a factor of 
~ 2 around the mid-point of reionization, highlighting that the 
scatter in the reionization redshifts of galaxies can be important 
in regulating star formation (e.g. Alva rez et al.|20 09 1. This homo- 
geneous reionization model further highlights the delay in the hy- 
drodynamic response of the gas to the UVB: M cr it only surpasses 
M coo i at z w 8, even though the ionizing background is effectively 
turned on at z IC = 9.3 (see also Oka moto et al.|2 008 1. 

It is natural to approximate this delay with the sound-crossing 
time-scale, t BC = 2_R v i r /c s (10 4 K) (which is comparable to the 
photo-evaporation time-scale, e.g. Haiman et al. 2001] | Shapiro | 
|et al.|2004) . For a halo of mass M, the sound-crossing time-scale 
can be approximated as (c.f. |Shapiro et al.(2 004): 

1/3 /-, , _\ -1 /n i.2\ - 1 / 3 

t sc « 200Myr | 



M 
1O 8 M 



1 + z 
10 



0.15 



(10) 

It would therefore be reasonable to expect M m in to evolve from 
M coo i at high redshifts to some Mo(T V i r ~ const) at low redshifts, 
with the transition occurring at z ~ z IC — Az sc [where Az sc is the 
redshift interval corresponding to eq. (10)], roughly over a time- 
scale scaling with the duration of reionization, Az rc . Indeed, we 
find that our results for M m j n are well fitted by: 



M min (z) = M coo i x 
with the transition function 

- 
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A/ coo l 



1 + exp 
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Az re 
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(12) 



We take Mq(z) to correspond to a fixed virial temperature To 
. ) x 10 4 K[]and define Az le to correspond to the redshift interval 
between Shi = 0.6 and Shi = 0.4. 



7 We caution that our simple ID simulations from Paper I are not well 
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Figure 3. Power spectrum A^ x of the ionization field at iehi = 0-75 (left), xjji = 0.50 (middle) and xjji = 0.25 (right) with different models of UVB 
feedback. Dot-Dashed: NF. Short-dashed: J=0.01. Dotted: .1=1. Dashed: IF. Long-dashed: SC. 
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Figure 4. Size distribution of the HII regions at iehi = 0.75 (left), ihi = 0.50 (middle) and rrjji 
Dot-Dashed: NF. Short-dashed: J=0.01. Dotted: J=l. Dashed: IF. Long-dashed: SC. 



In Fig. [6] we compare this analytic approximation (dotted lines) 
with the results from our simulation (dashed lines). We show the 
self-consistent model with the standard choice of the parameters 
(in this case the analytic formula and the result of simulation are 
indistinguishable). To test the robustness of the above approxima- 
tion, we also include two "extreme" models: (i) with molecular 
hydrogen cooled halos significantly contributing to reionization: 
Tcooi = 5 x 10 3 K (magenta curves); (ii) with a significantly higher 
emissivity £ = 100 (green curves). These runs correspond to z re = 
(9.3, 9.9, 12.0), Az re = (1.0, 1.3, 0.8) and Az sc = (2.0, 0.9, 2.5), for 
the fiducial, T coo i = 5 x 10 3 K, and £ = 100 models, respectively. 
Even in the extreme cases, our formula from eq. l | 1 I [ i provides a 
good fit. This general, physically-motivated expression facilitates 
implementing UVB feedback from an inhomogeneous reionization 
into models of galaxy formation. 



3.4 UVB Evolution 

In addition to the evolution of HII regions, UVB feedback affects 
the UVB itself. As star-formation is suppressed in vulnerable halos, 
the UVB decreases (with respect to the case without feedback). Our 
formalism explicitly computes the inhomogeneous UVB intensity, 
as described in j ]2.2.1| It therefore provides a useful prediction of 
the evolution of J21 during reionization. 

The UVB can be determined from measurements of the Lyman 



suited for low-redshift structure formation, where mergers and other 3D 
effects can be important. Hence, here we use a lower value of Mn than 
implied by extrapolating eq. |5J to z = 0; our choice is roughly consis- 
tent with cosmological simulations of homogeneous reionization jOkamoto] 
|et al.|2008| . 



0.25 (right) with different models of UVB feedback. 




Figure 5. Average value of M m i n versus redshift. Dot-dashed: no feed- 
back (NF). Short-dashed: our model with J21 = 0.01. Dotted: our model 
with J21 = 1. Dashed: instantaneous feedback (IF). Long-dashed self- 
consistent model (SC). Solid: for comparison we show M m j n (eq.[TJ as- 
suming a homogeneous reionization (HR) at zr\ = 9.3 (corresponding to 
the mid-point of reionization in the SC model). 



alpha forest in high-redshift quasar spectra. The ionization rate per 
baryon is found to b e remarkably constant from z ~ 2 — > 5, 
Tbkg ~ 10" 12 s" 1 ( jBolton & Haehnelt||2007| |Faucher-Giguere| 
|et al.|20 08 1. Then from z ~ 5 — > 6, there is evidence of a drop by a 
factor of few, though the uncertainty in the measurement is large at 
these high-redshifts when the Lya forest begins to saturate ([Bolton 



|& Haehnelt|2007l|Wyithe & Bolton|20TT] |Calverley et al.|201lf 

Such a drop could be caused by either an incomplete reionization 
(|Mesjnger 2010, McG reer et al.|20lT) , and/or by a rapid evolution 
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Figure 6. Average value of Af m ; n versus redshift. We compare the self- 
consistent model with the standard parameter choice (SC), with T coo i = 
5 X 10 3 K and with g = 100. Dashed: results of the simulation. Dotted: 
analytic expression in eq. |l 1| 
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Figure 7. Evolution of the UVB intensity J%\ (left-hand scale) and of 
the corresponding photoionization rate per baryon Tbkg (right-hand scale). 
Solid: average intensity { J)hii within HII regions in our SC model; the 
shaded region corresponds to the spread among HII regions. Dashed: aver- 
age intensity ( J)y in the entire simulation box. Dot-dashed: evolution of 
(-7) HII when feedback is neglected (NF). For comparison we show Tbkg 
at z ~ 6 (offset for clarity) as estimated from the Lyct forest ^Bolton &] 
|Haehnelt|2007]|Calverley et al.|2011) . 



in the abundance of LLSs ( |Furlanetto'& Oh 2005 , McQuinn et al. 
2011| >F1 



In Fig. [7] we show the redshift evolution of the UVB in our SC 
(solid curve) and NF (dot-dashed curve) models (^l/IYkg on the 
left/right vertical scales), together with z ~ 6 observational esti- 
mates by |Bolton & HaehneltH2007) and |Calverley et al.| ( pOTT| >. It 
is important to note that during UV-driven reionization, the UVB 
is expected to be bi-modal, with (J2i)ml inside HII regions and 
J21 ~ inside neutral regions, such that the volume-averaged 
UVB intensity is (J2i)v ~ (1 — xhi)(^2i)hii- This distinction 



8 Note that a large drop in the ionizing background must involve the photo- 
evaporation of LLSs (McQuinn et al. 201 1 1, since these systems regulate the 
progress of reionization during its final stages (e.g. Furlanetto & Mesinger 
|2009|[Croc"iani et al.|20lT]| Alvarez & Abel|2012) . We caution that our mod- 
els do not incorporate an evolving mean free path in ionized regions. There- 
fore, they likely underestimate the rise in the UVB in the late stages and 
following reionization. However, the relative impact of feedback is more 
robust to this uncertainty. 



is generally not made in Lyman forest studies, which thereby im- 
plicitly assume that reionization is over (though see Gallera m et al.| 
2006 and McGreer et al.|20iT|>, or that it is driven by X-rays with 
long mean free paths (e.g. Oh 2001 ). To highlight this point, in Fig. 
[Tjwe also show the evolution of (J2i)v in the SC model (dashed 
curve). 

From Fig. [7] we see that the average UVB intensity within HII 
regions in the SC model is quite flat, with a value of (J)hii — 
0.2, consistent with observational estimates at z ~ 6QComparing 
the solid and dot-dashed curves, we see that feedback in our SC 
model suppresses the UVB by a factor of < 2 towards the end of 
reionization. This factor is modest, compared to the uncertainties 
in the model parameters and observational measurements. 

The shaded region in Fig. [7] corresponds to the lcr scatter in 
(J2i)hii- This spread is roughly a factor of ~ 2 during most of 
reionization, decreasing to < 50% in the final overlap stages. Such 
a spread in the HH-region averaged intensity is roughly compara- 
ble to the overall spread in the intensity, as estimated by Mesinger 
|& Dijkstra| ( |2008^ and |Mesinger & Furlanetto| ( |2009^ who summed 
contributions from individual galaxies. This supports our formal- 
ism in §2.2. 1| and the corresponding simplification of using a 
roughly uniform ionizing background inside each HII region. 



3.5 Was Reionization "Self-Regulated"? 

Using the assumption of an instantaneous change in M m i n inside 
ionized regions (i.e. our IF model), many works claim that UVB 
feedback has a dramatic impact on reionization history and the cor- 
responding star-formation rate (e.g. Wyithe & Loeb|20 03 ; Onken 
& Miralda-Escud e|2004||Iliev et al.|2007||Wyithe & Morales|2007[ 
Wyithe & Cen|2007| |Kulkarni & Choudhury|2011[|Alvarez etal] 
2012| |Ahn et al.||2012| T Some works dub this effect to be "self- 
regulation", since when HII regions form, their subsequent growth 
is impaired. However, we show above that under more realistic pre- 
scriptions of UVB feedback, this effect is relatively minor (see also 
Mesinger & Dijkstra 2008 ). This qualitative claim is easy to under- 
stand: the time-scales associated with UVB feedback are compara- 
ble to the time-scale of reionization. Thus in the regimes we study, 
there is not enough time for UVB feedback to efficiently delay 
reionization. We now ask the question, "what are the physically- 
motivated regimes in which UVB feedback can in fact self-regulate 
reionization?" Naively one might expect that the strength of UVB 
feedback can be increased by (;) a higher UVB intensity; (ii) a 
more extended reionization; {Hi) a larger contribution of small ha- 
los, more susceptible to UVB feedback. We explore each of these 
in turn. 

In Fig. [8] we show the reionization history in models with dif- 
ferent ionizing photon efficiencies: £ =30, 50, 100 {left to right 
pairs of curves). We compare the evolution when UVB feedback 
is neglected (NF; dot-dashed line) with our self-consistent ap- 
proach (SC; long-dashed line). In fact, increasing the UVB in- 
tensity (J21 oc £), results in smaller feedback effects. This is 
because the dependence of M m i n on J21 is weak in our model 
(Af cr it oc J21 17 ), while for high values of £ , reionization proceeds 
faster. Therefore as the ionizing emissivity is increased, the de- 
crease in the duration of reionization is more relevant than the delay 



9 Choosing a lower/higher value of the ionizing efficiency, £, would de- 
lay/advance reionization and decrease/increase the resulting photoioniza- 
tion rate. The observational data still allow for a relatively broad range of 
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resulting from the increase in M cr it, and the evolution of reioniza- 
tion is driven by the emissivity. 

We next attempt to extend the duration of reionization, in order 
to give the UVB time to significantly suppress the gas content of 
galaxies. The most efficient way of significantly extending reion- 
izatiorp'lis to decrease the ionizing photon mean free path through 
the ionized IGM, R mip (e.g. |Alvarez & Abel|2012||Mesinger et al.| 
|2012[ >. Our fiducial model assumes J? m f p = 30 Mpc, consistent 
with observational and theoretical estimates at z « 6 (e.g.[Songaila 
|& Cowie|2010{|McQuinn et al.|20iT) . In Fig.|9]we show the reion- 
ization history of our fiducial model, but decreasing i? m f P to 5 Mpc. 
Although the reionization history is significantly extended, UVB 
feedback has a negligible effect, resulting in a delay of Az < 0.1. 
This is because the effective horizon imposed by 7? m f p restricts 
feedback to a very limited volume surrounding each galaxy. Reion- 
ization proceeds in a disjoint manner: with tiny, isolated HII regions 
emerging around newly-forming galaxies (Alvarez & Abel 2012 
|Mesinger et al.|2012) . Such a reionization scenario is not affected 
by UVB feedback. 

Finally, we turn to our point (iii): a larger contribution of small 
halos, more susceptible to UVB feedback. So far, we have consid- 
ered halos down to T coo i = 10 4 K: this is conservatively small, 
since winds and internal feedback can suppress star-formation in- 
side such halos (e.g. Springel & Hernquist 2003 , Pawli k & Schaye| 



2009 ; |Finlator et al.|201 1| >. However, some models argue for a sig- 
nificant contribution of galaxies hosted by even smaller halos (so- 
called minihalos), in which gas collapses through molecular cool- 
ing channels. In Fig. [5] we show the reionization history of such 
a toy model which allows halos down to T coo i = 5 x 10 3 K to 
emit ionizing photons with the same fiducial efficiency £ = 30. 
In this extreme case, UVB feedback has a dramatic impact on 
reionization history: delaying the advanced stages of reionization 
by Az ~ 1.5-2, and resulting in a slow evolution of xui during 
the first half of reionization. This result is qualitatively consistent 
with the recent results by |Ahn et al.| f2012) (who assumed instan- 
taneous UVB feedback, as well as instantaneous feedback from an 
H2-disassociative background; the later results in an even flatter, 
self-regulated regime during the early stages)! 1 " 1 "] 

We therefore conclude that reionization is self-regulating only if 
star-formation is efficient inside molecularly-cooled galaxies. We 
caution however that our implementation of UVB feedback from 
Paper I was calibrated only for atomic cooled halos. Furthermore, 
our models assume that the star-formation rate is proportional to 
the total amount of gas inside halos. We do not model the ability of 
gas to cool and form stars, which can be dramatically suppressed 
in minihalos by an H2-disassociative background before the bulk of 
reionization (assuming modest values for their ionizing emissivi- 
ties; e.g. Haiman et al.|1997| . 



10 Note that we already include star formation in halos down to the atomic- 
cooling threshold. The fractional abundance of these halos evolves more 
slowly than that of more massive halos, expected to have more efficient 
star-formation (e.g. Lidz et al. 2007). Therefore increasing the duration of 
reionization by changing the population of host halos is only possible if we 
include molecularly-cooled galaxies (see below). 

11 We caution that such simple models very likely overestimate the ion- 
izing efficiencies of minihalos, by ignoring the chemical and mechanical 
feedback from the first stars (e.g. Tornatore et al. 2007 , Whalen et al. 2008 1, 
X-ray pre-heating of the IGM (e.g. Ricotti & Ostriker 20041, as well as 
the baryon-DM velocity offset (e.g. Tseliakhovich & Hirata 2010). Simu- 
lations predict that these mechanisms are very efficient in suppressing star- 
formation in minihalos. 
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Figure 8. Evolution of the global neutral fraction xhi with different ioniz- 
ing efficiencies £ = 30, 50, 100. Dot-Dashed: no feedback. Long-dashed: 
self-consistent model. 
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Figure 9. Evolution of the global neutral fraction xm with B m f p = 
5 cMpc and with T v - lT = 5 X lO^ K. Dot-Dashed: no feedback. Long- 
dashed: self-consistent model. 

4 CONCLUSIONS 



Cosmic reionization inhomogeneously heated the IGM, affecting 
its cooling properties and photo-evaporating gas from the outskirts 
of galaxies. As a result, small dwarf galaxies inside the ionized 
patches would have a reduced gas reservoir for star-formation. This 
process involves a huge dynamic range; therefore large-scale reion- 
ization simulations generally either ignore this effect or include it 
with simple, fairly ad-hoc prescriptions. 

In this work, we use a tiered approach to study the impact of 
UVB feedback on reionization. Combining parameterized results 
for the minimum halo mass capable of hosting star-forming galax- 
ies (M m in, see Paper I) with semi-numerical simulations of reion- 
ization, we present large-scale (300 Mpc) reionization simulations 
which include UVB feedback in a physically-motivated manner. 
In our models, the ionizing emissivity of galaxies is assumed to 
be proportional to their baryonic content, and depends on the lo- 
cal values of the reionization redshift and the UVB intensity. These 
additions will be included in an upcoming version update of 2 1 CM- 
FAST. 

UVB feedback delays the end stages of reionization by Az < 
0.5, suppresses the large-scale ionization power by tens of percent 
(at fixed Shi), and decreases the mean photoionization rate by a 
factor of < 2 towards the end of reionization. These effects are 
quite modest, when compared to current astrophysical and obser- 
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vational uncertainties. By contrast, the popular assumption of an 
instantaneous change in M m i n results in a much larger effect, with 
a Az « 1.5 delay in reionization, and a factor of few higher aver- 
age M min . 

Assuming that UVB feedback acts on a time-scale roughly cor- 
responding to the halo sound-crossing time, we present an analytic 
formula for M m i n . This general expression fits our simulation re- 
sults very well, and can be easily included in models of galaxy 
formation. 

In our models, UVB feedback significantly delays reioniza- 
tion only if molecularly-cooled galaxies contribute significantly 
to reionization. Star-formation inside such galaxies is easily sup- 
pressed by mechanical feedback, X-ray pre-heating and the baryon- 
DM velocity offset. Hence, we conclude that UVB feedback is un- 
likely to self-regulate reionization. 

We thank Andrea Ferrara for comments on a draft version of this 
work. 
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